Anomaly in the heat capacity of Kondo superconductors 
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Using numerical renormalization group we study thermodynamic properties of a magnetic impurity described 
by the Anderson impurity model in a superconducting host material described by the BCS Hamiltonian. When 
the Kondo temperature in the normal state, Tk, is comparable to the critical temperature of the superconducting 
transition, T c , the magnetic doublet state may become degenerate with the Kondo singlet state, leading to a In 3 
peak in the temperature dependence of the impurity contribution to the entropy. This entropy increase translates 
into an anomalous feature in the heat capacity which might have already been experimentally observed. 
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I. INTRODUCTION 

After many decades of research, conventional supercon- 
ductors are still being actively studied both theoretically and 
experimentally 1 . An important experimental technique to 
study long-range-ordered states consists in determining what 
destroys that order 2 . Superconductivity can be suppressed, for 
example, by doping the host material with magnetic impuri- 
ties. Each magnetic impurity locally perturbs the supercon- 
ducting order; this can lead to a significant cumulative effect 
and ultimately, at large concentrations, to a complete suppres- 
sion of the superconductivity 2 *. 

In the normal state of a metal, the impurity magnetic mo- 
ment tends to be screened by the Kondo effect 3 -. In the su- 
perconducting state, however, the gap in the spectrum im- 
plies that there are no low-energy electrons to perform the 
screening 2 . The behavior thus depends on the ratio between 
the Kondo temperature Tk and the critical temperature T c for 
superconducting transition. In the Tjc <C T c limit, there is no 
Kondo effect since the gap is fully formed on the relevant en- 
ergy scale of Tk, while in the Tk 3> T c limit, Kondo physics 
does not affect the superconducting transition since the impu- 
rities are already fully nonmagnetic (Kondo screened) on the 
scale of T c . The non-trivial regime is thus Tjc ~ T c , when 
both effects compete. To study this competition and demon- 
strate the existence of the Kondo effect in superconductors, 
one can study, for example, the effects of the impurities on 
the density of states 4,5 or the thermodynamic 6 ' 7 and transport 
properties 8 . With advances in the tunneling spectroscopy, fo- 
cus has recently shifted to dynamic properties, i.e. variation 
of the local density of states induced by the presence of the 
impurity 2 i2iI9.. 

The effect of impurities on the host is the central focus in 
the field of quantum impurity physics^!. Studying thermo- 
dynamics of an impurity model is a common practice in this 
field, as one can extract many essential aspects of the im- 
purity behaviour (its magnetic properties, effective degrees 
of freedom, etc.). A reliable theoretical tool is the numeri- 
cal renormalisation group (NRG) 11 ' 12 ! 13 . The NRG has been 
applied, for example, to study the sub-gap bound states in 
isotropic 14 and anisotropic 15 Kondo model and in the Ander- 



son model 16 in both conventional and unconventional 17 *^ su- 
perconductors. Further NRG calculations focused on dynam- 
ical properties such as the spectral functio n 19 ' 20,21,22 and more 
recently on the transport properties of quantum dots attached 
to superconducting lead a 23,24,25 ' 26 ' 27 . Only few results of NRG 
calculations for thermodynamic quantities, such as the im- 
purity contribution to the magnetic susceptibility and the en- 
tropy, have been published so far^Si 2 !. In this paper we present 
a detailed investigation of these thermodynamic properties of 
an Anderson impurity in a BCS superconducting host. We 
find features in the thermodynamic behavior of magnetically 
doped superconductors (also known as the "Kondo supercon- 
ductors") in the dilute doping limit which appear to have been 
overlooked in theoretical treatments, but may have been al- 
ready experimentally observed. 



H. MODEL AND METHOD 

A conceptually proper way of studying the thermodynam- 
ics of a magnetic impurity in a superconducting host would 
start from a microscopic Hamiltonian for a conduction bath 
with explicit attractive electron-electron interaction (quartic) 
terms. The NRG can only be applied to non-interacting con- 
tinuum Hamiltonians or, more accurately, to Hamiltonians 
with at most quadratic (mean-field) terms. We thus focus on 
a mean-field BCS Hamiltonian, keeping in mind that this is 
only an effective low-temperature Hamiltonian and that its 
use in the NRG to compute temperature-dependent thermo- 
dynamic quantities is associated with a certain ambiguity. In 
a superconductor, the superconducting order is namely estab- 
lished gradually as the temperature is reduced below T c . In 
this work, we will make the following approximation: 



(0)(T) = (0(A = A(T)))(T), 



(1) 



i.e. at temperature T we take results computed in a calculation 
with a constant A chosen such that A = A(T). We will use a 
phenomenological approximation for the temperature depen- 
dence for the superconductor gap that is correct near T = 



2 



andT 



A(T) w 6 sc T c tanh 



(2) 



with <S SC = 1.76, a = 2/3, JC/CW = 1.43. 

Even for a single impurity, the order parameter becomes 
a spatially varying function A(r)^22. In this work, we dis- 
regard any spatial variation of the order parameter. In fact, 
we neglect any variation of A, considering it as a fixed (but 
doping-concentration dependent) quantity. In ppm to percent 
doping range, the distance between the impurities is much 
smaller than the coherence length in conventional supercon- 
ductors, thus the assumption of spatially constant order pa- 
rameter is justified when impurity-position-averaged thermo- 
dynamic properties are considered 30 . 

The impurity contribution to an expectation value is com- 
monly defined by taking the difference between the values for 
a doped and for a clean system: 



(0) imp = (O) - (O) 



(3) 



The impurity contribution can be further decomposed into two 
parts: the "intrinsic" local effect of an impurity and the aver- 
age effect of all the impurities on the bulk pairing order pa- 
rameter A (which needs to be calculated self-consistently). It 
is the first contribution that we are interested in in this work. 

We thus consider a single Anderson impurity embedded in 
a superconducting host described by the Hamiltonian H = 
Hhath + H\ mn + i/hyb- We model the fermionic bath with a 



'bath ~r -*-Mmp 

BCS model: 



#bath = 



E 



-feT 



(4) 



where Cka are conduction band operators, ek single-particle 
energies and A the BCS order parameter. For simplicity we 
assume that the band has a flat density of states in the nor- 
mal state, p(ui) — 1/2D, where 2D is the width of the 
conduction band. In the original BCS approach, the order 
parameter A has to be determined self-consistently through 
A = V 2~Z/c( c -fcl c /cT)' where V is some phenomenological 
effective coupling constant and the summation is to be cut off 
by the Debye energy u>d. In numerical calculations, however, 
we consider A to be some given parameter. The magnetic 
impurity Hamiltonian has the usual form 



H-, 



imp 



£](» 



i) 



f(n-l) 2 , 



(5) 



where e is the impurity energy level, U the electron-electron 
repulsion and n the level occupancy operator, n = dtd a - 
Finally, the hybridisation Hamiltonian is 



H hyh = vY,(dUka+H.c. 



(6) 



fc<7 



The coupling to the conduction band is characterized by the 
hybridisation strength T — TipV 2 . 



Due to strong exchange scattering of conduction band elec- 
trons on magnetic impurities 3 ^, magnetic impurity problems 
need to be solved using non-perturbative methods such as 
NRG. We performed the calculations with a discretization pa- 
rameter A = 3, four values of the twist parameter z and the 
truncation energy cutoff of -E C utoff = lOwjy, where lon is the 
characteristic energy-scale at the A-th NRG iteration. Only 
spin SU(2) symmetry can be explicitly used to simplify cal- 
culations with models of this class. 



III. NUMERICAL RESULTS 
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Figure 1: Thermodynamic properties of an Anderson impurity in a 
BCS superconductor, assuming a fixed value of A. The system is 
particle-hole symmetric: e = — [7/2. 

We present results for the following thermodynamic 
quantities: 1) the temperature-dependent impurity con- 
tribution to the magnetic susceptibility, Ximp(? 1 ) = 
{gfi B ) 2 /{k B T) ({S 2 Z ) - (Sf) ), where the subscript refers 
to the clean system (where the Hamiltonian is simply the band 
Hamiltonian iJband without the impurity terms), S z is the to- 
tal system spin, g the electronic gyromagnetic factor, ps the 
Bohr magneton and ks the Boltzmann's constant; 

2) the temperature-dependent impurity contribution to the 
entropy, S imp (T) = {E-F)/T-(E- F) /T, where E = 
(H) and F = -fc B TlnTr ( e - H ^ BT ); 

3) the temperature-dependent impurity contribution to 
the heat capacity, C imp (T) = k B ({H 2 ) - (H) 2 ) - 
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Figure 2: Close-up on the low-temperature region in the transition 
regime. 



k B «ff 2 ) - {Hf) Q . 

We first study the case of a constant gap parameter A, i.e. 
a fictitious gap parameter which does not change with the 
temperature. We computed thermodynamic properties for a 
range of A at fixed Kondo temperature, Fig.Q] (We use Wil- 
son's definition of T K , i.e. k B T K x(T K )/(gHB) 2 = 0.07 in 
the normal state.) For A <C Tk, we obtain the well-known 
results for an Anderson impurity in a normal metal exhibit- 
ing the formation of local moment at T ~ U and the Kondo 
screening at T ~ Tk, with the system ultimately ending up 
in the strong-coupling fixed point- . The presence of the gap 
is hardly reflected in the impurity thermodynamic properties, 
although we observe a small feature at T ~ A, when Kondo 
screening is interrupted and the remaining residual entropy 
is released. For A 3> Tk, the Kondo screening is again 
interrupted at T ~ A, but this time the system ends up in 
the local-moment fixed point rather than the strong-coupling 
fixed point: the impurity spin decouples entirely, so we obtain 
Curie-Weiss behaviour at low temperatures and a In 2 residual 
entropy. We bring attention to the small bump in the impu- 
rity entropy above In 2 at T ~ A, which is mirrored in an 
S-shaped feature in the impurity heat capacity. For A ~ Tk, 
the behavior is similar to that in the A 3> Tk regime, but the 
entropy increase at T ~ A becomes increasingly strong as 
the transition point is approached. In Fig. [2] we show a close- 
up on the transition region. Near the transition, the entropy 
rises to exactly In 3, forming a relatively broad plateau: at this 



point the magnetic doublet and the Kondo singlet states are 
degenerate, although this fixed point is unstable. The transi- 
tion is clearly first-order, since the quantum numbers of the 
ground state change. We present results only for the particle- 
hole symmetric case, as asymmetry does not affect the results 
qualitatively. 




Figure 3: Thermodynamic properties of an Anderson impurity in 
a BCS superconductor, taking into account the temperature depen- 
dence of the gap. 

We now turn to the case of BCS superconductors, where 
the gap is temperature dependent. The control quantity is in 
this case the ratio T c /Tk, which is adjusted in experiments 
by changing the magnetic impurity species and concentration 
(within the limits set by the miscibility of metals) or by chang- 
ing the doping x in the Kondo superconductor systems such 
as (Lai-zCe^A^. 

As previously described, we compute the thermodynamic 
quantity (O) (T) at temperature T from the NRG results for a 
fixed value of gap A = A(T). The results of such a calcula- 
tion for a range of the Tk /T c ratios are shown in Fig. [3] For 
T > T c , the results are equivalent to those for a magnetic im- 
purity in a normal metal. For T < T c , the gap starts to open. 
On a logarithmic scale the features in thermodynamic quan- 
tities appear sharp (cusp-like), but a close-up reveals that the 
change of slope is in fact continuous, albeit rapid. The results 
are qualitatively equivalent to those obtained for fixed A; a 
notable difference, however, is the absence of the In 3 plateau. 
The plateau is replaced by a rounded peak, since the condition 
of exact degeneracy of the three many-particle states is satis- 
fied only at a single temperature. The peak in the entropy is 
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then reflected in a characteristic temperature variation of the 
impurity heat capacity, which has a positive and negative part: 
this is the characteristic signature to be sought after in exper- 
iments. We emphasize that the impurity heat capacity may be 
negative, since it is defined as the difference of two total heat 
capacities (which themselves are, of course, positive quanti- 
ties). 

IV. DISCUSSION AND CONCLUSION 

We find that in Kondo superconductor systems an anomaly 
appears in the impurity heat capacity for Tk ~ T c . This is 
a large effect, since the overall variation in the impurity heat 
capacity is of the order of ks and it occurs in a temperature 
interval of approximately one decade of the temperature scale. 
The positive part of the anomaly occurs at low temperatures, 
i.e. in the exponential tail of the BCS heat capacity, thus the 
heat-capacity enhancement is sizable and should easily be ob- 
servable already for impurity concentrations even below one 
permil. The negative part of the anomaly occurs in the tem- 
perature range where the heat capacity of a superconductor 
is enhanced, thus the heat-capacity reduction would be more 
difficult to measure. Nevertheless, the magnitude of the ef- 
fect (0.1 — 0.2/cb per impurity) is large enough that this effect 
could be detected at somewhat higher impurity concentration 
in the permil range. 

At high doping levels, the simple picture of indepen- 
dent impurities is no longer valid due to interactions be- 
tween the impurities. Owing to their pair breaking proper- 



ties, magnetic impurities lead to a suppression of the super- 
conducting transition temperature T ^ 2 ' 32 ' 33 ' 34 ' 35 . In both the 
Abrikosov-Gorkov treatment 33 with classical impurity spins 
and the Nagaoka-Suhl-approximation treatment by Zittarz and 
Miiller-Hartmann 4 ' 35 with quantum impurity spins, the de- 
crease of the critical temperature is found to be proportional 
to the concentration divided by fcsp(0), where p(0) is the 
density of states at the Fermi level in a pure normal-state 
metal. For impurity concentrations in the permil range, ele- 
mental superconductors with relatively high T c such as lead or 
niobium are still superconducting. Another class of systems 
where these effects may be detected are Kondo superconduc- 
tor alloys. In fact, excess specific heat at low temperatures 
had already been observed in (Lai_ x Ce x )Al2 alloy long age 6 . 
Given the progress in both experimental and theoretical tech- 
niques, it would be interesting to revisit this problem focusing 
on the low-doping regime to detect the predicted reduction of 
the heat capacity before the onset of the enhancement and to 
look for possible semi-quantitative agreement. We conclude 
by observing that a magnetic impurity in a bath may be con- 
sidered as an open quantum system thus this problem is also of 
interest in the context of thermodynamics of nanosystems 36 . 
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